APPLICATIONS OF CONVEX ANALYSIS TO THE 
SMALLEST INTERSECTING BALL PROBLEM 
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Abstract: The smallest enclosing circle problem asks for the circle of smallest radius enclosing a 
given set of finite points on the plane. This problem was introduced in the 19th century by Sylvester 
[17]. After more than a century, the problem remains very active. This paper is the continuation 
of our effort in shedding new light to classical geometry problems using advanced tools of convex 
analysis and optimization. We propose and study the following generalized version of the smallest 
enclosing circle problem: given a finite number of nonempty closed convex sets in a reflexive Banach 
space, find a ball with the smallest radius that intersects all of the sets. 
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1 Introduction and Problem Formulation 

A more general form of the smallest enclosing circle problem is called the smallest enclosing 
ball problem. Given a set P = {pi,...,p n }, n > 1, on a Banach space X, it is always 
possible to find a ball IB (a; r) such that 

P C B(a;r). 

The smallest enclosing ball problem asks for such a ball with the smallest radius. 
Consider the following optimization problem 

minimize f(x), x £ X, (1-1) 

where the function / therein is defined by 

f(x) = max{| \x — pi 1 1 : i = 1, . . . , n}. 

The smallest enclosing ball can be found by solving (1.1). 

Numerous articles have been written to study the smallest enclosing ball problem as 
well as its generalizations from both numerical and theoretical viewpoints. The reader are 
referred to [6, 7, 16, 18] and the reference therein for recent developments as well as the 
history of the problem. 
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In this paper we propose and study a new 
problem called the smallest intersecting ball 
problem as follows: given a finite number 
of nonempty closed convex sets in a Banach 
space, find a ball with the smallest radius 
that intersects all of the sets. It is obvious 
that when the sets under consideration are 
singleton, the problem reduces to the small- 
est enclosing ball problem. This is the contin- 
uation of our previous work from [15, 14] in 
an effort to shed new light to classical geom- 
etry problems using new tools of nonsmooth 
analysis. 

Let fij, i = 1, . . . ,n, n > 1, be nonempty 
closed convex sets in a Banach space X. Let 




Figure 1: A Smallest Intersecting Ball Problem. 

x be any point in X. Then there always exists r > such that 



B(x; r) n fi; ^ for alii = 1, . . . , n. 
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We are looking for a ball with the smallest radius r > (if exists) such that property (1.2) 
holds. Define the function 



D(x) = m&x{d(x; fij) : i = 1, . . . , n}, 
where the distance function generated by a set Q is given by 

d(x; O) = inf { | |a; — cj\\ : oj € O}. 



;i.3) 



(1.4) 



It is not hard to see that the function D is convex since each function d{x; fij), i = 1, . . . , n, 
is convex. Consider the problem of minimizing the function D on X below 



minimize D(x) subject to x £ X. 



(1.5) 



As we will see in Section 3, if X is a reflexive space and Oj, i = 1, . . . ,n, have no point in 
common, such a smallest intersecting ball can be found by solving problem (1.5). For this 
reason, we are going to use the following standing assumptions throughout the paper unless 
otherwise specified: 

X is a reflexive Banach space and Oj, i = 1, . . . , n, n > \, are nonempty closed convex sets 
in X with 

n? =1 fii = 0. 

Notice that if the Oj, i = 1, . . . ,n, have a common point, then the smallest intersecting 
ball problem has no solution (unless balls of radius are allowed) . 

The smallest intersecting ball problem is an example of facility location problems. In 
contrast to most of the existing facility location problems which deal with locations of 
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negligible sizes (points), this new problem deals with those that involve locations of non- 
negligible sizes (sets). The problem is obviously mathematically interesting with promising 
applications to location models in which it is possible to access the entire of each location 
from a point in it. The difficulty when dealing with the smallest intersecting ball problem 
comes from the nonsmooth nature of the cost function D in (1.5), especially when the norm 
in A is Non-Euclidean. Our approach in this paper is to study the problem from both 
theoretical and numerical viewpoints using new tools from convex analysis and optimization. 
The results we are going to present in this paper provide improvements and generalizations 
of many results in [6, 7]. 

We organize the paper as follows. Section 2 provides necessary tools from convex analysis 
and optimization for solving the smallest intersecting ball problem. In Section 3, we study 
the problem from theoretical aspects. Section 4 is devoted to developing an algorithm 
of subgradient type to solve the smallest intersecting problem in finite dimensions. The 
MATLAB implementations of the algorithm are also presented. 

2 Tools of Convex Analysis 

This section provides important constructions and results from convex analysis that will be 
used in the next sections. Most of the material presented here can be found in [2, 9, 19]. 

Let A be a normed space with the dual space A*. A function / : A — >• M is called 
convex if 

/(Ax + (1 - X)y) < A/(x) + (1 - X)f(y) for all x, y G A and A € (0, 1). (2.1) 

If the inequality in (2.1) becomes strict for x ^ y, the function / is called strictly convex. 
A subset O of A is called convex if 

Ax + (1 - X)y € U for all x, y G and A € (0, 1). 

It is not hard to prove that a nonempty closed subset O of A is convex if and only if the 
corresponding distance function (1.4) is convex. Note that the distance function f(x) = 
d(x; f2) is Lipschitz continuous on A with modulus one, i.e., 

\f{x) - f(y)\ <\\x -y\\ for all x, y G A. 

Let (•, •) be the dual pair between A and A*. An element x* G A* is called a subgradient 
of a convex function / at x if the following holds 

(x*, x — x) < f(x) — f(x) for all x G A. 

The set of all subgradients of / at x is called the subdifferential of / at x denoted by df{x). 

Convex functions and subdifferentials have several important properties as far as opti- 
mization is concerned. For instance, a convex function / has a local minimum at x if and 
only if it has an absolute minimum at x. Furthermore, the following generalized version of 
the Fermat rule holds 

x is a minimizer of / if and only if G df(x). (2-2) 
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It is well known that the subdifferential of the distance function (1.4) has a close connection 
to the normal cone of the generating set fi. Recall that the normal cone of a convex set 
at x € SI is defined by 

N(x;Sl) = {v G X* : (v,x-x)<0 for all x € fi}. (2.3) 

The projection from a point x 6 X to a set fi is 

n(i;(]) = {w£!l : = d(x;fi)}. (2.4) 

The following representation of subdifferential for the distance function (1.4) will play an 
important role in the subsequent sections of the paper. The proof of the formulas can be 
found in [5], while their various extensions are presented in [13]. 

Proposition 2.1 Let Q be a nonempty closed convex set of a Banach space X and let 
x € X. Suppose that n(cc;fi) ^ (which is always the case when X is reflexive). Then 



dd(x; fi) 



dp(x -u)n N(u>; fi) if x <£ Q, 
N(x;Q)nB* ifxeSl, 



where IB* is the closed unit ball of X* and lo is any element ofIl(x;Q). 
In particular, if X is a Hilbert space, then U(x;£l) is singleton and 



dd(x; O) 



N(x;Q)nlB if xeQ, 



In this proposition we also observe that when x € 0, one has Tl(x;£l) = {x}. Since 
dp(0) = IB*, the formula 

dd(x; Q) = dp(x -u)n N{Cu; Q) 

holds for any x £ X. 

Finally, we present the following well-known subdifferential rule that involves "max" 
functions. 

Theorem 2.2 Let X be a Banach space and let fi : X — > M, i = 1, ... ,n, be continuous 
convex functions. Define 

f(x) = max{/i(x) :i = l,...,n}. 

Then 

8 f(x)= co {dfi(x):i£l(x)}, 
where I(x) = {i = 1, . . . , n : f(x) = fi(x)}. 
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3 The Smallest Intersecting Ball Problem: Theoretical 
Aspects 

This section is devoted to theoretical analysis of the smallest intersecting ball problem. 
We are able to provide improvements and generalizations of many results in [6, 7]. Our 
approach is based mostly on tools of convex analysis and optimization. 

The following proposition allows us to reduce the smallest intersecting ball problem to 
a nonsmooth convex optimization problem in the reflexive space setting. For this reason, 
we will identify the smallest intersecting ball problem with problem (1.5). 

Proposition 3.1 Consider the minimization problem (1.5). Then x 6 X is an optimal 
solution of this problem with r = D{x) if and only if H3(x;r) is a smallest ball that satisfies 
(1.2). 

Proof: Suppose that x is an optimal solution of (1.5) with r = D{x). Since Qi, i = 1, . . . , n, 
have no point in common as in the standing assumptions, one has 

D{x) = mi{D(x) : x G X] = r > 0. 

This implies 

d(x; Qi) < r for all i = 1, . . . , n. 
Since X is reflexive, there exist u)j € fij, i = 1, . . . , n, satisfying 

\\% — < r. 

It follows that Qi 6 JB(x; r) PI Oj, and hence 

B(x; r) n Qi ^ for alH = 1, . . . , n. 

Suppose there exists r' < r and x! 6 X with 

B(x'] r') n Cli ^ for alH = 1, . . . , n. 

Then 

d(x ; Q,j) < r < r for all i = 1, . . . , n. 

This implies D(x') < r 1 < r = D(x), which is a contradiction. Thus S3(x;r) is a smallest 
ball we are looking for. 

We are now going to justify the converse. Let us first prove that r = D{x). Since 
]B(x; r) (~l Qi ^ 0, one has 

d(x; Qi) < r for all i = 1, . . . , n. 

This implies D(x) < r. Assume by contradiction that D(x) < r. Let r' satisfy D(x) < r' < 
r. Then 

B(x; r') n Qi ± for all i = 1, . . . , n. 
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This contradicts the minimal property of r. Thus r = D{x). Let x be any point in X and 
let r' = D(x). Then 

B(x; r') n 0^ + for all i = 1, . . . , n. 
This implies r < r' or D{x) < D{x'). Therefore, x is an optimal solution of (1.5). □ 

In what follows, we will prove that under natural assumptions on the sets 0^, i = 1, . . . , n, 
such a smallest intersecting ball does exist. We are going to use the fact that on a Banach 
space, any convex lower semicontinuous function is weakly lower semicontinuous. 

Proposition 3.2 Suppose that there exists i = 1, ... ,n, such that fij is bounded. Then the 
smallest intersecting ball problem (1.5) has a solution. 

Proof: Without loss of generality, suppose that £l\ is bounded. Define 

r = inf{D(x) : x G X}. 
Let (x n ) be a minimizing sequence for problem (1.5). That means 

D(x n ) — > r as n — > oo. 

Let N G IN with 

d(x n ; fi x ) < D(x n ) < r + 1 for all n > N. 
Then there exists a sequence (ui n ) in such that 

\\x n — L) n \\ < r + 1 for all n > N. 

Since (uj n ) is a bounded sequence, (x n ) is also bounded. As X is reflexive, there exists a 
subsequence (x„J that converges weakly to x. This implies 

D(x) < liminf D(x nk ) < r 

because D is weakly lower semicontinuous. Therefore, x is a solution of problem (1.5). □ 

Proposition 3.2 implies that the smallest enclosing ball problem (1.1) always has a 
solution because each Jlj = {wj}, i = 1, . . . ,n, is obviously bounded. However, in general, 
the smallest intersecting ball problem (1.5) may not have any solution. 

Example 3.3 Let X = M 2 with the Euclidean norm. Consider Oi = {0} x M and 

^2 = {(x,y) G M 2 : y > -,x > 0}. 

x 

Then the smallest intersecting ball problem (1.5) generated by and does not have 
any solution. 

In the case where the smallest intersecting ball (1.5) has a solution, the solution may 
not be unique as shown in the example below. 

Example 3.4 Let X = M 2 with the Euclidean norm. Consider Q\ = {(x, y) G M 2 : y > 1} 
and fii = {(x,y) £ M 2 : y < -1}. Then any x G M x {0} is a solution of the smallest 
intersecting ball problem (1.5) generated by S7i and ^2- 
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Example 3.5 Consider X = M 2 with the "max" norm ||(xi,X2)|| = max{|xi|, |rc2 1 } - Then 
the ball H3(x; r) in X, where x = (x\,X2) and r > 0, is the square 

S(x; r) = [x\ — r, x\ + r] x [x2 — r, x% + r]. 

Problem (1.5) can be equivalently interpreted as follow: find a smallest square S(x;r) that 
intersects fij for alH = 1, . . . , n. Using different norms on X, we obtain different intersecting 
ball problems. 

Lemma 3.6 Let X be a Hilbert space and let uJi G X , i = 1, . . . ,n. Then the function 

s(x) = max{||x — cjj|| 2 : i = 1, . . . ,n}, n > 1, 

is strictly convex. 

Proof: We are going to prove that for x ^ y and t £ (0,1), one has 

s(tx + (1 - t)y) < ts{x) + (1 - t)s{y). 

By induction, we only need to show that the function p(x) = \ \x\\ 2 is strictly convex and the 
function g(x) = m&x{gi(x), g2(x)} is strictly convex if both g\ and 52 are strictly convex 
functions. Indeed, for t € (0, 1) and x,y S X, one has 

p(tx + (1 - t)y) = \\tx + (1 - t)y\\ 2 

= t 2 \\x\\ 2 + 2i(l - t)(x,y) + (1 - t) 2 ||y|| 2 
<t 2 ||x|| 2 + 2t(l-t)||x||.||y|| + (l-i) 2 ||y|| 2 
<t 2 ||x|| 2 + t(l-t)(||x|| 2 + |y|| 2 ) + (l-t) 2 ||y|| 2 
= t||x|| 2 + (1 - t)\\y\\ 2 = tp{x) + (1 - t)p(y). 

Notice that the equality holds if and only if ||x|| = ||y|| and (x,y) = ||x||||y||. This implies 
ll x — y\\ 2 = an d hence x = y. Therefore, p is strictly convex. 
Now let t G (0, 1) and x ^ y. Then 

9l {tx + (1 - t)y) < t 9i (x) + (1 - t) 9i (y) < tg{x) + (1 - t)g{y) for i = 1, 2. 

This implies 

g (tx + (l - t )y) < tg(x) + (1 - t)g{y). 

The proof is now complete. □ 

The following proposition gives an example of a smallest intersecting ball problem which 
has a unique solution. We will use a natural convention that JB(c; 0) = {c} for any c G X. 

Proposition 3.7 Let X be a Hilbert space and let r > 0. Suppose that fi, = JB{oji\r), 
i = 1, . . . , n, are closed balls in a Hilbert space X. Then the smallest intersecting ball problem 
(1.5) generated by Qi, i = l,...,n, has a unique solution. Moreover, this unique solution 
coincides with the unique solution of the smallest enclosing ball problem (1.1) generated by 
the centers of the balls {uji}, i = 1, . . . ,n. 
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Proof: Let us first show that in this case the function D in (1.3) has the following repre- 
sentation 

D{x) = max{||x — cJi\\ : i = 1, . . . ,n} — r. (3-1) 

Indeed, let 

J(x) = {i G 1, . . . , n : x ^ fij}. 

Since fij, i = 1, . . . ,n, have no point in common by the standing assumptions, J(x) ^ 0. 
For any i S {!,... J(x) and for any j £ J(x), one has 



|x — Wj| | < r < \\x — u)j\ 



It follows that 



D{x) = max{d(x; Q,) : i = 1, . . . , n} 
= max{d(x; fij) : i G «/(^)} 
= max{||:c — cjjll : i £ J(x)} — r 
= max{||j; — Wi|| : i = 1, . . . ,n} — r. 

Thus (3.1) has been justified. Using representation (3.1), we see that x is a solution of 
problem (1.5) if and only if it is a solution of the minimization problem 

minimize s(x) = max{||x — Wj|| 2 : i = 1, . . . , n}, x € X. (3.2) 

Since s is strictly convex by Lemma 3.6, problem (3.2) has a unique solution. Therefore, 
problem (1.5) also has a unique solution. Notice that a; is a solution of problem (3.2) if 
and only if it is the solution of the smallest enclosing ball problem (1.1) generated by {wj}, 
i = 1, . . . , n. The proof is now complete. □ 

Example 3.8 In M 2 with the Euclidean norm, consider the balls f2i = iB((0, 3);3), = 
1B((-2,0);1), and 3 = {B((2, 0); 1)}. Then x = (0,0) is the solution of problem (1.5) 
generated by these balls, but this solution is the solution of the smallest enclosing ball 
problem (1.1) generated by the centers of the balls. 

In what follows we are going to prove that in the Hilbert space setting, the smallest 
intersecting ball problem (1.5) generated by closed balls with different radii also has a 
unique solution although the solution may not coincide with the solution of the smallest 
enclosing ball problem generated by their centers. 

Proposition 3.9 Let X be a Hilbert space and let £li = JB(ui;ri), ri > 0, i = 1, . . . , n, be 
closed balls in a Hilbert space X. Then the smallest intersecting ball problem (1-5) generated 
by Q,i, i = 1, . . . , n, has a unique solution. 

Proof: Similar to the proof of Proposition 3.7, one has 

D(x) = maxjUx — u>i\ \ — r^ : i = 1, . . . , n} for all x € X. 
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Define 

Pi(x) = \\x - -n,i= 1, . . . ,n. 

Let 

£ = — max{rj : i = 1, . . . , n}. 
Then pi(x) > £ for all i = 1, . . . , n and for all x £ X. Consider the optimization problem 

minimize h(x) = max{(j>j(x) — I) 2 : i = 1, . . . , n}. (3-3) 

Notice that Pi(x) — £ > for all i = 1, . . . , n and for all x G X. Then it is not hard to see 
that x is a solution of the smallest intersecting ball problem (1.5) if and only if it is also a 
solution of problem (3.3). Similar to the proof of Proposition 3.7, one sees that the function 
h in (3.3) is strictly convex and hence problem (3.3) has a unique solution. The proof is 
complete. □ 

For each x G X, the set of active indices for D at x is defined by 
I(x) = {ie{l,...,n}: D(x) = d(x; J^)}, 

and let 

Ai(x) = dp(x - uii) n N(uji; Qi). 

where oji £ n(x;Oj). Notice that the definition of Ai(x) does not depend on the choice of 
Ui by Proposition 2.1. It is also clear from the definition that I(x) ^ for any x G X. 
Moreover, if i G I(x), then D(x) > because n" =1 f2j = as in the standing assumptions, 
and hence 

d(x;Sli) = D(x) > 0. 

This implies x ^ 

Proposition 3.10 Consider the smallest intersecting ball problem (1.5). Then x G X is 
an optimal solution of the problem if and only if 

G co {Ai(x) : i G I{x)}. 

Proof: It follows from Theorem 2.2 that 

dD(x) = co {dd(x; : % G I(x)} 
= co {Ai(x) : i G I(x)}. 

The result then follows from the subdifferential Fermat rule (2.2). □ 

Corollary 3.11 Let X be a Hilbert space. Consider the smallest intersecting ball problem 
(1.5). Then x is a solution of the problem if and only if 

x G co{wj : i G I(x)}, (3.4) 

where uji = H(x; Qi). 

In particular, if Qi = {aj}, i = 1, . . . ,n, then x is the solution of the smallest enclosing 
ball problem (1.1) generated by a{, i = 1, . . . ,n, if and only if 

x G co{aj : i G -f(x)}. (3-5) 
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Proof: According to Proposition 3.10, the element x is a solution of the smallest intersecting 
ball problem (1.5) if and only if 

G co {Ai(x) : i G I(x)}. 
For each i G I(x), one has x £ f2j, and hence 

Ai(i) = ^i^)^ = ^^f^ 

It follows that 

G co{Aj(x) : i G I(x)} 
if and only if there exists Aj > 0, i G /(cc), such that XneiYs) ^« = 1 an d 

o= y 

^ D(x) 

This equation is equivalent to 

0= ^2 Xi(x - Qi) or x = ^ ^i^i, 
i<=I(x) i£l(x) 

which is equivalent to (3.4). 

Notice that (3.4) is equivalent to (3.5) when Qj = {a^}, % = 1, . . . ,n. The proof is now 
complete. □ 

We say that the smallest ball JB(x; r) touches a target set f2j, % = 1, . . . , n, if fijin JB{x; r) 
is singleton. 

Corollary 3.12 Let X be a Hilbert space. Consider the smallest intersecting ball problem 
(1-5). Then any smallest intersecting ball touches at least two sets among f2j, i = 1, . . . ,n. 

Proof: Let us first prove that \I(x)\ > 2 if x is a solution of problem (1.5). Suppose by 
contradiction that I(x) = {io}- Then by (3.4), 

x G Qi and d(x; O^ ) = D(x) > 0. 

This is a contradiction. Let us now show if i G I(x), then JB(x;r) touches fij, where 
r = D{x). Indeed, in this case d(x; ilj) = r. If there are u,v G such that 

u,v £ U3{x; r) D Q^, u ^ v. 

Then 

1 1 ti — 1 1 < r = d(x; Jlj) and ||u — x|| < r = d(x; fij). 

It follows that u, u G II(x;f2j). This is a contradiction because II(x;f2j) is singleton. Thus 
EB(x;r) touches Oj. The proof is now complete. □. 

It is obvious that co {<2j : i G L(x)} C co {cij : i = 1, . . . , n}. Thus our result in Corollary 
(3.11) covers [7, Theorem 3.6]. It is also possible to prove that the solution of the smallest 
intersecting ball problem (1.5) generated by closed balls in a Hilbert space belongs to the 
convex hull of their centers as in the next proposition. 
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Proposition 3.13 Let X be a Hilbert space. Suppose that = JB(uf,ri), ri > 0, i = 
1, . . . ,n, are closed balls in X. Let x be the unique solution of problem (1-5). Then 

x € co {uii : i G L(x)}. 

Proof: Let x be the solution of the smallest intersecting ball problem (1.5) generated by 
fij, i = l,...,n. Since x for all i € L(x), by (3.4) from Corollary 3.11, there exist 

Aj > 0, i 6 L(x), such that Yliei(x) = I and 



EAi(wi + ri-^z — -\). 



X — 

iG/(i) 11 



For any i € I(x), one has 1 1 as — = r + rj, where r is the radius of the smallest intersecting 
ball. It follows that 

This implies 

x = t > — uji G co {uii : i G 

v-v Xi .^r + n 11 v >s 

The proof is now complete. □ 

Example 3.14 Let aj,i = 1,2,3, be three points in M 2 with the Euclidean norm and let 
x, be the solution of problem (1.1) generated by a^, i = 1,2,3. By Corollary 3.12, one 
has |/(2;) | = 2 or |/(x)| = 3. If |/(x)| = 2, say L(x) = {2,3}, then x, € co{a2,a3} and 
||x — 0.2 1 1 = \\x — «3 1 1 by Corollary 3.11. In this case 

«2 + a 3 



This also implies (02 — a±, a$ — a\) < 0. Conversely, if (02 — 01, 03 — 02) < 0, then the angle 
of the triangle formed by a\, 02, and 03 at vertex a\ is obtuse (we allow the case where a^, 

i = 1,2, 3, are on a straight line). One can easily see that J( Q2 — — -) = {2, 3}, and ° 2 Q3 

02 + O3 

satisfies the assumption of Corollary 3.11. Then x = because of the uniqueness of 

the solution. In this case we have \L{x)\ = 2. 

Thus |/(x)| = 2 if and only if one of the angles of the triangle formed by 01,02, and 03 
is obtuse. In this case the solution of problem (1.5) is the midpoint of the side opposite to 
the obtuse vertex. 

If none of the angles of the triangle formed by 01,02,03 is obtuse, then \L(x)\ = 3. In 
this case, x is the unique point that satisfies 

X £ CO {fill, 02, 03}, \ \x — Ol| I = \\x — 02! I = \ \x — 03! I, 

or x is the center of the circumscribing circle of the triangle. 
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Figure 2: Smallest Intersecting Ball Problems for Three Disks in M 2 



Let us now consider the solution of problem (1.5) with the target sets being three disjoint 
disks in M 2 . 

Example 3.15 Let Oj = lB(cji,ri), r\ > 0, i = 1,2,3, be disjoint disks in ]R 2 with the 
Euclidean norm. We use bd(£li) to denote the boundary of Qi, which is the circle of center 
u>i and radius n, i = 1, 2, 3. Let x be the unique solution of the problem. 

Let us consider the first case where one of the line segments connecting two of the centers 
intersects the other disks. For instance, the line segment connecting 0J2 and uj 3 intersects 
fix. Let U2 = <^2^3 H bd(Q2) and u 3 = U2OJ3 n bd(Qs). Let x be the midpoint of U2U3. Then 
I(x) = {2, 3} and we can apply Corollary 3.11 to see that x is the solution of the problem. 

Now we only need to consider the case where any line segment connecting two centers 
of the disks does not intersect the remaining disk. Let 

u\ = u\uj2 n bd(Qi), v\ = uJTcJs n bd(Qi), 

u 2 = U2U3 n bd(tt 2 ), V2 = n bd(fl 2 ), 

u 3 = TJ3UJ1 n bd(fl s ), v 3 = u 3 u) 2 n bd(Q 3 ) 

a\ = uTprrTi n bd(Q\), 02 = ^2 m 2 H ^2, 03 = ^3^3 H ^3, 

where m\ is the midpoint of U2V3, 1^2 1S the midpoint of u 3 vi, and 7713 is the midpoint of 
u±V2- If one of the angles: v^a~\v 3 , u^a^vi, U\a^V2 is greater than or equal to 90°. For 
instance, if U2a\v 3 is greater than or equal to 90°. Then J(mi) = {2,3}, and x = m\ is 
the unique solution of the problem by Corollary 3.11. Now if all of the afore- mentioned 
angles are acute, then I{x) = 3 and the smallest disk we are looking for is the unique disk 
that touches three other disks. The construction of this disk is the celebrated problem of 
Apollonius; see, e.g., [8]. 

We are going prove that a smallest intersecting ball generated by n convex sets, n > 1, 
in a M m can be determined by at most m + 1 sets among them; see Figure 1 for the 
visualization of this property. The proof is based on a known results for points (see [18, 
Lemma 1 (iii)]), which can be easily proved by Corollary 3.5 and the Caratheodory theorem 
[10, Corollary 1, Sec. 3.5]. 
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Lemma 3.16 Let P = {p\, . . . ,p n }, n > 1, be a set of finite points in M m with the Eu- 
clidean norm and let lB(x;r) be the smallest enclosing ball for problem (1.1) generated by 
points in P. Then there exists a subset Q C P and 2 < \Q\ < m + 1 such that lB(x;r) is 
also the smallest enclosing ball of problem (1.1) generated by points in Q. 

Proof: By Corollary (3.11), one has 

I £ co {pi : i G I(x)}. 

By the Caratheodory theorem, there exists an index set J C I(x) with \J\ < m + 1 and 

x G co {pj : j G J}. (3.6) 

It is clear that | J\ > 2 because n > 1. Let Q = {pj : j G J}. Then 2 < \Q\ < m + 1 
and \\x — q\ \ = r for all q € Q. By converse of Corollary 3.11, one has that !B(x;r) is the 
smallest enclosing ball of problem (1.1) generated by points in Q. The proof is complete. 
□ 

Proposition 3.17 Let X = M m with the Euclidean norm. Consider problem (1.5) in 
which Qi, i = l,...,n, are disjoint. Suppose that lB(x;r) is a smallest intersecting ball of 
the problem. Then there exists an index set J with 2 < |Jj < m + 1 such that lB(x;r) is 
also a smallest intersecting ball of problem (1.5) in which the target sets are ftj, j € J. 

Proof: Let lB(x;r) be a smallest intersecting ball of problem (1.5) with target sets fij, 
i = 1, ... ,n. By Corollary 3.11 and Corollary 3.12, one has 

x € co {iOi : i G L(x)}, 

where \L(x)\ > 2 and G 11(5;; Qi). Again, by Corollary 3.11, !B(x\ r) is the solution of the 
smallest enclosing ball (1.1) generated by {cuj : i G L(x)}. 

Applying Lemma 3.16, one finds a subset J C L(x), 2 < \ J\ < m + 1 such that lB(x; r) 
is the solution of the smallest enclosing ball generated by {ojj : j G J}. Then 

x G co {uji : i G J} 

and \\x — 0Jj\\ = r for all j G J. Now consider problem (1.5) generated by {IX,- : j G J}. 
Applying Corollary 3.11, we also see that x is a solution of this problem because 

x G co {ujj : j G J}, 

where ojj G H(x;Qj) and d{x;VLj) = r for each j G J. Moreover, JB(x;r) is a smallest 
intersecting ball for the problem. The proof is now complete. □ 

Now we are going to give a generalization of [7, Theorem 4.4]. Our approach, which is 
based on the proof of [10, Proposition 1, Sec. 10.2], allows us to give an estimate of the 
radius of the smallest intersecting ball for problem (1.5) generated by closed balls in M m . 
Notice that [7, Theorem 4.4] holds in M 2 for the classical smallest enclosing circle problem. 
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Theorem 3.18 Let X = M m with the Euclidean norm. Consider the smallest intersecting 
ball problem (1.5) generated by the closed balls = JB(u)f,ri), i = 1,... ,n. Let 

r min = min{ri : i = 1, . . . , n}, r max := max{rj : i = 1, . . . , n}, 
£ = min{m + 1, n}, P = {uji : i = 1, . . . , n}, 

and Zei -B* = JB(x;r) be the smallest enclosing ball. Then 



In particular, 



2(£ - 1) 

c/iam P—max {diamQi : i = 1, . . . , n} < diam -B* < y diam P—min {diamQi : i = 1, 

Proof: For any i, j = 1, . . . ,n, one has 

ll^i — < ||x — + ||x — ojj\\ < 2r + n + < 2r + 2r ma:r . 

Thus 

diam P < 2r + 2r max 

and the first inequality in (3.7) holds true. 

Let us prove the second inequality of (3.7). By the Caratheodory theorem [10, Corol- 
lary 1, Sec. 3.5] and the proof of Proposition 3.13, there exist k < min{m + l,I(x)} < 
min{m + 1, n} and Aj > 0, % = 1, . . . , k, Yli=i \ = 1 (we reorder the indices if necessary) 
such that 

fiiUJi, where m = ■ — , fj, = } — ■ 

i=i i=i 

and 

| \x — 0Ji\\ = r + n,i = 1, . . . , k. 

One has 

— 0Jj\\ 2 = \\0Ji — x\\ 2 + \ \ujj — x\\ 2 — 2{bJi — X,UJj — x) 

= (r + r-j) 2 + (r + rj) 2 - 2(u)i - x,ujj - x). 

Thus 

k k k 

- ujjW 2 = (r + rj) 2 + } j fJ,j(r + r,) 2 - 2(^J/^u;j - x,Uj - x) 

i=l i=l i=l 

k 

= (r + r-j) 2 + ^ l JL ii r + r i) 2 - 
i=i 

It follows that 

k k 



V>iHj\\ui -Uj\\ 2 = 2^/t 8 (r + nY 

i,j=l i=l 
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We also have 

k k k 



This implies 



Ms = (2^^) - 2^ ^ 1 - ^ = nr - - ~T' 

ij=l,i^j i=l i=l 



k k £ — I 

^ fJ>iHj\\ui - Wj|| 2 = 2^/ii(r + rj) 2 < (diam P) 2 — —. 

i,j=l i=l 



Vi(r + ri) 2 = - ^ (r + rj) 2 = — (r + r*), 
fin+r p, 



k 

^(r + r,) <^(diam P) 2 
i=l " 

Using the formula for /i, we arrive at 

1 1 — 1 

(r + r min ) 2 < -— — (diam P) 2 . 
This implies (3.7). The second estimates follow from (3.7). The proof is complete. □ 



Since 



one has 



4 Subgradient Algorithm and Its Implementation 

In this section let (X, || • ||) be a normed space where X = M m and let p(x) = \ \x\ \ be the 
norm function on X. We are going to present and justify an algorithm of subgradient type 
to solve problem (1.5) numerically and illustrate its implementations using MATLAB. 

Theorem 4.1 Let i = 1, ... ,n, be nonempty closed convex subsets of X such that at 
least one of them is bounded. Picking a sequence {oik} of positive numbers and a starting 
point x\ € X, consider the iterative algorithm: 

x k+ i = x k - a k x* k , k E IN. (4.1) 

Let the vectors x* k in (4-1) be given by 

x* k e dp(x k - uj k ) n N{uj k ; a), (4.2) 

where uj k € n(x k ; fij) and i is any index chosen from the following index set 

H x k) = {i = l,...,n: D(x k ) = d(x k ; fij)}. 

Assume that the given sequence {a k } in (4.1) satisfies the conditions 

oo oo 

^2a k = oo and ^^a|<oo. (4.3) 

k=l k=l 
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Then the iterative sequence {x k } in (4.2) converges to an optimal solution of the smallest 
intersecting ball problem (1.5) and the value sequence 

V k = mm{D( X] ):j = l,...,k} (4.4) 

converges to the optimal value V in this problem. 
Furthermore, we have the estimate 

where £ 2 = X^fcLi a L an< ^ d-2,(x\;S) denotes the distance generated by the Euclidean norm 
from x\ to the solution set S of the problem. 

Proof: By Proposition 3.2 the smallest intersecting ball problem under consideration has a 
solution. Observe that the function D in (1.5) satisfies a Lipschitz condition with Lipschitz 
constant k = 1. We have 

dD(x k ) = co{dd(x k ;Qi) : i £ I{x k )} 

= co {dp(x k - u k ) (~1 N(u k ; Oj) : i £ I(x k )}, 

where uj k € Tl(x k ;£li). Notice that under the standing assumptions, x k ^ for i € I(x k ). 
It follows that for any i 6 I(x k ) one has 

9p(a; fe - Wfc) n iV(wfe; Oj) C dD(x k ). 

We also have that dd(x k ; f2j) = dp(x k — uj k ) n N(uj k ; fij) is nonempty. Since all norms in X 
are equivalent, it suffices to show that 

I \x k — x\ I2 — > and 14 — > V, 

where || • H2 is the Euclidean norm in X and x is a solution of problem (1.5). However, these 
follow directly from the well-known results on the subgradient method for convex functions 
in the so-called "square summable but not summable case"; see, e.g., [4, 1]. □ 

One important features of the subgradient method is that the subgradient x^. for each 
k is not uniquely defined. This also reflects in the following direct consequence of Theorem 
4.1. 

Corollary 4.2 Let X = M m with the Euclidean norm. Consider the smallest intersecting 
ball problem (1.5). For each k E JN, the subgradient x* k in Theorem 4-1 is computed by 

* _ x k - u k 

%h 11 111 



\x k - U) k \ 



where uj k = H(x k ; fij) and i is an index chosen from I{x k ). 

In particular, if£l% = JB(a;ri), i = 1, . . . ,n, are closed balls in X. Then the subgradient 
x* k has the following explicit representation 

x* k = -p^ for an index i E I{x k ). 

I \ x k c% 1 1 
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-6 -4 -2 2 4 6 8 10 



MATLAB RESULT 


k 




T T 

v k 


1 


(2,2) 


10.29563 


1,000 


(-1.05477,3.05560) 


7.13408 


10,000 


(-1.05562,3.05553) 


7.13408 


100,000 


(-1.05555,3.05555) 


7.13408 


400,000 


(-1.05556,3.05556) 


7.13408 


800,000 


(-1.05556,3.05556) 


7.13408 


1,000,000 


(-1.05556,3.05556) 


7.13408 



Figure 3: A Smallest Intersecting Ball with Euclidean Norm to Square Targets. 

Example 4.3 Consider X = M 2 with the Euclidean norm. The target sets are the squares 
fij = S(uji;ri), u)i = (wii,(jj2i), i = l,...,n, where 

S(oJi',ri) = [u) U ~ ri,u) U + x [u 2i - ri,u 2i + r»]. 

Let the vertices of the ith square be denoted by vu = (uJu + ri,uj2i + ri),V2i = (uju — ri,uj2i + 
n),V3i = {uiu - ri,uj2i - n),V4i = (cJii + ri,u)2i - n) and let Xk = (xik,x 2 k)- Fix an index 
i € I(xk). Then the vectors in Theorem 4.1 are given by 

if xu - > Ti and x 2 k - w 2 ; > n, 
if xik - wu < -n and x 2 k - u 2l > n, 
if xik - u)u < —n and x 2 k - u 2l < -r,, 
if xik - wi, > n and x 2k - u 2 i < ~n, 



x k 


- v u 


\\Xk 


-vu\\ 


x k 


- V 2 i 


\\Xk 


~V2i\\ 


Xk 


- v 3i 


\\xk 


-VSi\\ 


Xk 


- V 4 i 



Xi, = < 



\\Xk- Vju\ 

(0,1) 

(o,-i) 

(1,0) 
(-1,0) 



if |arifc -U)u\ < fi and x 2k - Lu 2i > n, 
if \xi k -L0u\ < Ti and x 2k - Lu 2i < -r< 
if xik - wu > n and \x 2k —u> 2 i\ < n, 
if x\k - wii < -n and \x 2k - uj 2 A < 



It is also not hard to determine the index set I(xk) sequence and Vk for each k E IV. Thus 
the algorithm is explicit. 

Consider the target sets 0», i = 1, . . . , 7, to be the squares with centers (—8, 8), (—7, 0), 
(-4,-1), (2,0), (2,-6), (7,1), and (6,5) and the radii n = {1,2,3,0.5,2,1,1} for i = 
1, . . . , 7, respectively. A MATLAB program is performed for the sequence a& = 1/k satis- 
fying (4.3) and the starting point x\ = (2,2); see Figure 3. 

Observe that the numerical results computed up to five decimal places yield an optimal 
solution x ~ (—1.05556,3.05556) and the optimal value V ~ 7.13408. 
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When working with a norm in X that is different from the Euclidean norm, it may be 
difficult to find the distance functions, the projections to sets, as well as the subdifferential 
of the norm. The following remark allows us to have an intuitive way to find a subgradient 
xt, k E JV, in Theorem (4.1) in the case X = M 2 with the "sum" norm. 

Remark 4.4 Let X = M 2 with the "sum" norm p(x) = \x±\ + \x2\, x = (xi,X2). The ball 
]B(x;t), x = (xi,X2), t > 0, is the following diamond shape 

]B(x; t) = {(xi, x 2 ) € X : \xi - x\\ + \x 2 - x 2 | < t}. 

The distance from x to a nonempty closed set Q and the corresponding projection are given 
by 

d(x;fi) = min{t > : B(x;t) Hfl ^ 0} (4.5) 

and 

Il(x; O) = iB(x; i) n ft, where t = d(x; O). 
Moreover, the subdifferential dp(x), i£l, has the following explicit representation 



9p(xi,x 2 ) = 



' [-1,1] x [-1,1], 


if 


(x u x 2 ) = (0,0), 


[-1.1] X{1}, 


if 


X\ = 0, X2 > 0, 


[-l,l]x{-l}, 


if 


X\ = 0, X,2 < 0, 


{l}x[-l,l], 


if 


X\ > 0, X2 = 0, 


{-l}x [-1,1], 


if 


X\ < 0, §2 = 0, 


{1}x{1}, 


if 


£x > 0, 82 > 0, 


{I}x{-1}, 


if 


X\ > 0, X2 < 0, 


{-1}x{1}, 


if 


5l < 0, X2 > 0, 


{-I}x{-1}, 


if 


5l < 0, X2 < 0. 



By considering the "sum" norm in X, we are able to introduce a new smallest intersecting 
ball problem in which a "ball" is a diamond shape. The algorithm is going to be implemented 
in the following example. 

Example 4.5 Let us consider an example when X = M 2 with the "sum" norm. Let f2j 
be the squares 5(cjj;rj), i = l,...,n, given in Example 4.3. Notice that a ball ]B(x;r), 
x = (xi,X2), in X is the diamond shape 

IB(x; r) = {(xi, X2) 6 M 2 : \xi — xi| + \x2 — x~i\ < r}. 

Therefore, the smallest intersecting ball problem (1.5) can be interpreted as follows: find a 
diamond shape in M 2 that intersects all n given squares. Using the same notation for the 
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vertices of the target set f2j as in Example 4.3, one can see that the vectors x£ in Theorem 
4.1 are given by 



Xi, = < 



' (1,1), 


if 


Xlk - 


UJli 


> n and X2k - 


ui2i > n, 


(-1,1), 


if 


Xlk ~ 


UJli 


< —ri and x 2 k 


- uj 2i > r i} 


(-1,-1), 


if 


Xlk - 


Wli 


< — r 4 and x 2 k 


- uj 2i < -n 


(1,-1), 


if 


Xlk ~ 


UJli 


> r l and x 2k ~ 


u) 2 i < -n, 


(0,1), 


if 


\xik - 


- U)t 


< n and X2k - 


- uj 2i > n, 


(0,-1), 


if 


\xik - 


- u± 


< r 4 and x 2 k - 


- ui 2i < -r i} 


(1,0), 


if 


Xlk ~ 


UJli 


> ^ and \x 2 k - 


- ui 2l \ < n, 


. (-1,0), 


if 


Xlk - 


UJli 


< and \x 2k 


- uj 2i \ < ri. 



Consider the target sets = 1,... ,7, to be the squares with centers (—5,3), (—3,0), 
(—2, —3), (0, —8), (4, —3), (3, 0), and (5, 4) and the radii n = 1 for i = 1, . . . , 6, respectively. 
A MATLAB program is performed for the sequence = 1/k satisfying (4.3) and the 
starting point x\ = (2,0); see Figure 4. 




-4 -2 2 4 



MATLAB RESULT 


k 


Xk 


v k 


1 


(2,0) 


8 


10 


(0.67500,-0.17897) 


6.85397 


100 


(0.49999,-0.25853) 


6.75040 


1,000 


(0.50000,-0.25049) 


6.75040 


10,000 


(0.50000,-0.25004) 


6.75004 


100,000 


(0.50000,-0.25000) 


6.75000 


200,000 


(0.50000,-0.25000) 


6.75000 



Figure 4: A Smallest Intersecting Ball with Sum Norm to Square Targets. 

Observe that the numerical results computed up to five decimal places yield optimal 
solution x w (0.50000, -0.25000) and the optimal value V « 6.75000. 

Similar observations for X = M 2 with the "max" norm can be easily seen: 

Remark 4.6 Let X = M 2 with the "max" norm p(x) = max{|xi|, |x2|}, x = (x±,X2)- The 
ball S3(x;t), x = (x±,X2), t > 0, is the following square 

lB(x; t) = [xi — r, x\ + r] x [x 2 — r, X2 + r]. 

The distance from x to a nonempty closed set Q and the corresponding projection are given 
by 

d(x; 0) = mm{t > : JB(x; t)f]Q^$} 
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II(x; CI) = B(x; t) n Cl, where f = d(x; Cl). 
Moreover, the subdifferential dp(x), x G X, has the following explicit representation 



dp(xi,x 2 ) = < 





N - 




<i} 






if 


(xi,x a ) = (0,0) 


{(0,1)} 












if 


|*i | < x 2 , 


{(o,-i)} 












if 


x 2 < -\xi\, 


{(1,0)} 












if 


X\ > \x 2 \, 


{(-i,o)} 












if 


x\ < -\x 2 \, 


{(vi, v 2 ) G M 2 


M- 


i-N 


= 1, vi > 0, 


v 2 


>0} 


if 


x\ — x 2 > 0, 


{(«1,«2) G JR 2 


M- 


i-N 


= 1, vi > 0, 


v 2 


<0} 


if 


x\ = —x 2 > 0, 


{(«i,w 2 ) G fi 2 


M- 


t-M 


= 1, vi < 0, 


r 2 


<0} 


if 


xi = x 2 < 0, 


{(«1,W 2 ) G JR 2 


K'i - 


i-N 


= 1, vi < 0, 


l'2 


>0} 


if 


X\ = — .T 2 < 0. 



Example 4.7 Let us consider an example when X = M 2 with the "max" norm. Let Oj 
be the squares given in Example 4.3. Notice that a ball S3(x;r), x = (xi,X2), in X is the 
square 

lB(x; r) = [x\ — r, x\ + r] x [x2 — r, xi + r]. 

Therefore, the smallest intersecting ball problem (1.5) can be interpreted as follows: find a 
smallest square in M 2 that intersects all n given squares. Using the same notation for the 
vertices of the target set Cli as Example 4.3, one can see that the vectors x* k in Theorem 4.1 
are given by 



(1,0), 


if 


\x 2 k 


- U 2 i 


< 


Xlk 


- Uli 


and 


Xlk > Wli 


+ r 


(-1,0), 


if 


\x 2 k 


- U 2 i 


< 




- Xlk 


and 


Xlk < Wit 


- r 


(0,1), 


if 


\xik 


- 


< 


Xlk 


- u 2i 


and 


X 2 k > UJ 2i 


+ r 


(0,-1), 


if 


\xik 


- w Xi 


< 


u 2l 


- x 2k 


and 


X 2 k < u 2i 


- r 



where i € I(xk). The sequence (Vfc) is determined based on D(xk)- Fix any i G I{xk). 
Then 





(wii + r ? ), 


if 


\x 2 k 


- 


< 


Xli 


- Wu, 


Xlk > uu - 


f n, 


(wit - 


- r 4 ) - xife, 


if 


\x 2 k 


- u 2i 


< 




- Xlfc) 


Xlk < 


- r», 


x 2k — 


(w 2l +r»), 


if 


\xik 


- wu 


< 


X 2 k 


- w 2l , 


X 2 k > u 2l 


+ r», 




- n) - a; 2 fc, 


if 


\Xlk 


- Wii 


< 


u 2 t 


- x 2k , 


x 2 k < u2i 


- n 
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Consider the target sets = 1,...,6, to be the squares with centers (—5,7), (—2,0), 
(2,-5), (7,-2), (3,2), and (7,8) and the radii n = {1,1,0.5,1,2,0.5} for i = 1,...,6, 
respectively. A MATLAB program is performed for the sequence a& = 1/k satisfying 
(4.3) and the starting point x\ = (-2,3); see Figure 5. Observe that the numerical results 




MATLAB RESULT 


k 


Xk 


v k 


1 


(-2,3) 


8.50000 


10 


(-0.52500,1.54603) 


7.04603 


100 


(0.02973,1.00990) 


6.50000 


1,000 


(0.02973,1.00100) 


6.50000 


10,000 


(0.02973,1.00010) 


6.50000 


100,000 


(0.02973,1.00000) 


6.50000 


200,000 


(0.02973,1.00000) 


6.50000 



Figure 5: A Smallest Intersecting Ball with Max Norm to Square Targets. 

computed up to five decimal places yield an optimal solution x ~ (0.02973, 1.00000) and 
the optimal value V ~ 6.50000. 

The advantage of the algorithm comes from the fact that we are able to deal with the 
smallest intersecting ball problem generated by target sets of different types and different 
norms. Although a faster subgradient algorithm may be applied to this problem, we have 
chosen the simplest one for demonstrations. 
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